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i 1 Abstract 

X 

To help constrain the algorithms used in reconstructing high-energy muon events incident on the Sudbury Neutrino Observatory 
1 (SNO), a muon tracking system was installed. The system consisted of four planes of wire chambers, which were triggered by 
O scintillator panels. The system was integrated with SNO's main data acquisition system and took data for a total of 95 live days. 
^ Using cosmic-ray events reconstructed in both the wire chambers and in SNO's water Cherenkov detector, the external muon 
i " i tracking system was able to constrain the uncertainty on the muon direction to better than 0.6°. 



^ 1. Introduction 
(N 
<N 

{N| The Sudbury Neutrino Observatory (SNO) was a large wa- 
,— I ter Cherenkov detector optimized for detecting solar neutri- 
l/") nos created from the 8 B reaction in the main pp fusion chain. 

In addition to solar neutrinos, the Sudbury Neutrino Observa- 
,— ' tory was also sensitive to high-energy muons that traverse the 
I volume of the detector. A small fraction of these events are 
£> neutrino-induced muons from atmospheric neutrinos, while the 
k> large remaining fraction come from cosmic rays created in the 
j_j upper atmosphere. It is possible to discriminate between these 
C$ muon sources by looking at the angular distribution of incom- 
ing muons. The combination of large depth and the relatively 
fiat topography in the vicinity of the detector attenuates al- 
most all cosmic ray muons entering the detector at zenith angle 
cos (6 Z ) > 0.4. The study of muon events in the SNO detec- 
tor provides measurements of the absolute flux of atmospheric 
neutrinos and constraints on the atmospheric neutrino mixing 
parameters Am| 3 and 023 IFI- While the latter measurement is 
more strongly constrained by other experiments [2, 3], the for- 
mer is unique to the SNO experiment. 
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To facilitate a clean measurement of the zenith distribution 
of muons entering the SNO fiducial volume, an accurate un- 
derstanding of the muon reconstruction algorithm is necessary. 
This includes both the angular and spatial resolution of high- 
energy muons which enter the detector. Determining the accu- 
racy of the muon tracking reconstruction algorithm, however, 
relies almost entirely on Monte Carlo simulations. Although 
the detector response to muons was benchmarked against se- 
lected cosmic -ray data, there is not an external calibration 
source that can provide a consistency-check to the accuracy of 
the reconstruction algorithm. This is in sharp contrast to the 
case for SNO's response to neutrons and low energy electrons, 
which was calibrated with multiple sources to a precision of 

~ i% m. 

We present in this paper a means by which the SNO experi- 
ment was able to calibrate its muon tracking algorithm via the 
use of an external muon tracking system. The External Muon 
System (EMuS) allowed SNO to simultaneously reconstruct se- 
lected cosmic-ray events in two independent systems, thereby 
providing a cross-check on the tracking algorithm. The EMuS 
experiment ran for a total of 94.6 live days during the last phase 
of the SNO experiment. 

This paper is divided as follows: Section|2]describes the main 
SNO experiment, Section [3] describes the SNO muon recon- 
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struction algorithm, Section HI describes the characteristics of 
the EMuS apparatus, Sectionp] describes the criterion for ac- 
cepting events, and finally Section[6]discusses the analysis used 
to calibrate the SNO tracking algorithm against data taken with 
the EMuS system. 

2. The Sudbury Neutrino Observatory 

The SNO detector consisted of a 12-meter-diameter acrylic 
sphere filled with 1 kiloton of D2O. The 5.5-cm-thick acrylic 
vessel was surrounded by 7.4 kilotons of ultra-pure H2O en- 
cased within a barrel-shaped cavity, 34 m in height and 22 
m in diameter. A 17.8-meter-diameter geodesic structure sur- 
rounded the acrylic vessel and supported 9456 20-cm-diameter 
photomultiplier tubes (PMTs) pointed toward the center of the 
detector. A non-imaging light concentrator was mounted on 
each PMT to increase the effective photocathode coverage to 
54% [5|. The detector is described in detail elsewhere [6|. 

SNO was located in the Vale Creighton mine in Ontario, 
Canada at a depth of 2.092 km (5890 ± 94 meters water equiva- 
lent) with a flat overburden. At this depth, the muon rate in- 
cident over the geodesic sphere and integrated over the sea- 
sonal variation is 62.9 + 0.2 /i/day across an impact area of 216 
m 2 JTl. Muons entering the detector produce Cherenkov light at 
an angle of 42° with respect to the propagation direction of the 
muon. Cherenkov light and light from delta rays produced by 
the muon illuminate an average of 5500 PMTs, whose charge 
and timing information are recorded. The amplitude and timing 
response of the PMTs were calibrated in situ using a light dif- 
fusing sphere illuminated by a laser at six distinct wavelengths 
J3). This laser ball calibration was of particular relevance to 
the muon fitter because it provides a timing and charge calibra- 
tion for multiple photon hits on a single PMT. Other calibration 
sources used in SNO are described elsewhere J6]|7l. 

Data taking in the SNO experiment was subdivided into three 
distinct phases for measurement of the solar neutrino flux. In 
the first phase, the experiment ran with pure D2O only. The 
solar neutral current reaction was observed by detecting the 
6.25 MeV y-ray following the capture of the neutron by a 
deuteron. For the second phase of data taking, approximately 
0.2% by weight of purified NaCl was added to the D2O vol- 
ume to enhance the sensitivity to neutrons via their capture 
on 35 C1. In the third and final phase of the experiment, 40 
discrete 3 He and 4 He proportional tubes were inserted within 
the fiducial volume of the detector. This enhanced the neutron 
capture cross-section to make an independent measurement of 
the neutron flux, by observing neutron capture on 3 He in the 
proportional counters. Results from the measurements of the 
solar neutrino flux for these phases have been reported else- 
where H M EB H3 ECS EEU • 



3. Muon Reconstruction with the SNO Detector 

The SNO muon reconstruction algorithm fits for a through- 
going muon track based on the charge, timing, and spatial dis- 
tribution of triggered PMTs. Using a maximum likelihood 



3 0.03 

C 



025 



0.02 



0.015 



0.01 



0.005 



7 photoelectrons 



J 




20 photoelectrons 



& ''Mi; 



20 25 30 
Charge (GSC) 

Figure 1: The normalized PMT charge distribution measured in (scaled) 
pedestal-subtracted ADC charge for the case of 7 and 20 photoelectrons striking 
a single PMT. The red line indicates the prediction from the charge parameter- 
ization model used in the reconstruction. 



method, the fitter is able to determine a variety of muon tracking 
parameters, including the muon's propagation direction, impact 
parameter with respect to the center of SNO, the total deposited 
energy, and a timing offset. The likelihood is defined as: 



PMTs 



L= f[ J]PNWdPQ(Qi\n)PT(m 



n=\ 



(1) 



where n is the number of detected photons, P N (n\Aj) is the prob- 
ability of n photoelectrons being detected for A; expected num- 
ber of detected photoelectrons, Pg(2/|") is the probability of 
seeing charge Qi given n photon hits, and iV(f,-|n) is the proba- 
bility of observing a PMT trigger at time t given n photon hits. 

The heart of the fitter lies in the first probability term, which 
is calculated based on Monte Carlo simulations. Muons were 
simulated at discrete impact parameter values with random di- 
rections through the detector. These simulations were used to 
create lookup tables for how many photoelectrons are expected 
to be detected by a PMT at a given position with respect to a 
muon track with a given impact parameter. 

The second term further refines the fit by including the charge 
information from the PMTs, and allows an estimate of the total 
energy deposited by the muon, correcting for offline PMTs and 
the neck of the detector. This probability was calculated by sim- 
ulating multiple photon hits on all of the PMTs in SNO. For a 
given number of photon hits, the resulting charge distribution is 
modeled as an asymmetric Gaussian with the widths extracted 
from simulations. This fit model agrees well with the simula- 
tions for many photon hits, and acceptably for few photon hits 
(see Figure[T|. 

The third term in the likelihood refines the fit by including the 
PMT timing. For each PMT, the time residual can be calculated 
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Figure 2: The angular difference (as defined in Eq. [5J of Monte-Carlo muon 
tracks through the SNO detector (solid histogram). The angular distribution is 
fit to the function outlined in Eq.[4](solid line). The results from the fit are given 
in TableQ] 



as: 



tres - tpMT.i ~ to - 



CD 



(2) 



where tpMTj is the recorded time on a given PMT, to is the time 
offset term in the likelihood fit, d\ is the distance the muon 
travels within the detector before emitting the Cherenkov pho- 
ton, c is the speed of light in vacuum, di is the distance the 
Cherenkov photon traveled, and cd is the average speed of light 
in D2O/H2O medium (21.8 cm/ns). The Cherenkov photon is 
assumed to have an angle of 42° with respect to the muon track, 
making d\ and c/2 well-defined. The probability of the time 
residual is modeled as a Gaussian centered at zero with cor- 
rections to include estimates of prepulsing and late light as a 
function of the number of photon hits. 

The SNO muon fitter maximizes the likelihood function for 
the impact parameter, direction, deposited energy, and timing 
offset using the method of simulated annealing with downhill 
simplex lfl6l . After determining the parameters that maximize 
the likelihood, a set of data quality measurements are used for 
background rejection. 

The muon fitter is found to have good reconstruction accu- 
racy on simulated muons. Figure [2] shows the angle (8 mr ) be- 
tween the Monte Carlo generated muon direction (u g ) and the 
reconstructed muon direction (it r ): 



@mr — COS (Wg • Wy) 

This is fit to a weighted double Gaussian function: 



p(0) = AO 



fe-^+(l-f)e- 



(3) 



(4) 



Figure 3: The impact parameter difference of Monte-Carlo muon tracks 
through the SNO detector (solid histogram). The distribution is fit to the func- 
tion outlined in Eq.|5](solid line). The results from the fit are given in Table[T] 
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Table 1: Accuracy of the muon fitter based on Monte Carlo simulations. Fit 
parameters for mean (//), widths (cr and mcr), and relative weight (1 — f) are 
given in Equations|4]and|5] 



The additional f?-dependence is introduced in order to ac- 
count for the phase space available. 

The fit parameters are summarized in Table [T] Although the 
tails are non-Gaussian, this fit gives a reasonable estimate for 
the uncertainty for the angular resolution. Figure [3] shows the 
impact parameter reconstruction accuracy. The distribution is 
fit to the sum of two Gaussians: 



p(x) = A 



fe 



+ (!-/)* 



(5) 



with the fit parameters also summarized in Table [T] Monte 
Carlo studies show that the reconstruction accuracy of the muon 
direction and impact parameter are uncorrected. 

4. The External Muon System 

The External Muon System consists of a series of 128 single- 
wire chambers arranged into four planes and triggered by three 
large scintillator panels (see Figure|4]i. The wire chamber cells 
and electronics were provided by the University of Indiana. 
Each cell is 7.5 cm wide and has a square cross-section with 
the corners trimmed into a near-octagonal shape. The cells are 
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Figure 4: Diagram of the EMuS detector. See the text for more details. 
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Figure 5: Diagram of the EMuS electronics system. See the text for more 
details. 



2.564 m in length and possess a single 50 //m diameter tungsten 
wire running through the center. The wire is held at a positive 
potential of 2500 V (2700 V) while running on the surface (un- 
derground) for electron drift and collection. A gas mixture of 
90%Ar-10%CO2 was used in order to achieve high efficiency 
and stability, and to meet safety regulations for underground 
operations. 

When a muon passes through the system, it deposits energy 
in the scintillator and ionizes atoms in each of the wire cham- 
bers it passes through. The scintillator converts the energy into 
light that is then detected by PMTs in a fast process (~ns). In 
the wire chambers, the high voltage draws the ionization elec- 
trons to the wire in a slow drift process (~ //s). The drift time 
is proportional to the closest distance between the muon track 
and the wire, allowing track reconstruction using timing and 
position. The measured drift time for each wire is the time dif- 
ference between when the scintillator fired and when the drift 
electrons reached the wire. 

The scintillator consists of three large rectangular panels 
(350 x 70 x 5 cm 3 ) which cover the active region of the EMuS 
detector. The panels were acquired from the KARMEN neu- 
trino experiment IfTTl . and consisted of Bicron BC412 scintilla- 
tor read out at each end by four Photonis XP2262 PMTs. The 
signals from the PMTs were sent to a LeCroy 2249A Analog to 
Digital Converter (ADC) and a discriminator. If both ends of a 
panel fire in coincidence, a start signal was sent to the wire read- 
out modules, and the ADC modules recorded the pulse-height 
of eachPMT. 

Each wire chamber was monitored by an individual Front- 
End Electronics (FEE) card which outputted an ECL signal if 
a pulse is detected on the wire. The ECL signal was sent to a 
LeCroy 3377 Time to Digital Converter (TDC) with a readout 
window of 4. 1 [is. In order to mitigate high levels of electronic 
noise in the pre-amplifiers, the readout cables were sent through 
an additional ECL-NIM-ECL converter (see Figure[5]l. 

The EMuS system was deployed on the deck of the SNO 
experiment, 12 m above and 3 m west of the center of the de- 
tector. Due to space and solid-angle considerations, the planes 
were inclined at a 55° from horizontal. A survey was per- 
formed to determine the position of each of the wires with re- 
spect to the SNO detector. The dominant sources of uncertainty 
associated with the wire positions relative to the SNO detec- 
tor are summarized in Table [2] The largest uncertainty stems 
from determining X-Y coordinates of the EMuS detector. By 
comparing survey results with other known location markers at 
the detector, the X-Y coordinate was determined to better than 
+0.53 cm. The reference point used for the Z-coordinate of the 
detector was only known to +0.32 cm, and thus added as an 
uncertainty to the EMuS location. Other uncertainties on the 
locations of the wires included uncertainties on the floor level, 
the placement of the wires within the modules, the spacing be- 
tween wires, and the gaps between the modules. These addi- 
tional uncertainties do not apply equally to all wires, and have 
a maximum combined value of +0.30 cm. The final uncertainty 
on the SNO-EMuS coordinate translation based on this survey 
was ±0.68 cm. 
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Figure 6: The drift time for simulated electrons inside the EMuS wire chambers 
plotted as a function of starting radius. The plot shows drift times for both cir- 
cular (boxes) and octagonal (circles) cross-sectional geometries. The quadratic 
fit (solid line) is accurate to within 5% at the maximum simulated radius. 



4.1. Time to Radius Conversion 




R (cm) 

Figure 7: Drift time as a function of radius for data taken at Bates Laboratory 
(surface measurement). The color axis indicates the number of events that re- 
construct with the given radius and time. The vertical error bars are Garfield 
simulations of the drift time. 



Well-determined models of electron drift and diffusion in a 
gas 1 18 1 predict that the timing of a wire chamber hit with re- 
spect to the scintillator trigger can be used to measure the dis- 
tance of closest approach of the muon. This time-to-radius con- 
version function, r(f), has been simulated and measured for the 
EMuS system. 

The Garfield gas simulation IT9] was used to generate ex- 
pected r(t) curves as a function of gas pressure and applied 
voltage. The code was not able to perfectly model the shape 
of the wire chambers so two similar geometries were used to 
check the effects of this imperfect modeling: a circle with ra- 
dius 3.75 cm, and a regular octagon with a longest radius of 
4.06 cm. Simulated electrons were generated at 10 points along 
the longest radius, and the mean drift time for each point was 
calculated. Figure [6] shows that the two r(t) curves agree to 
within 2%. A parabolic fit to this data is accurate to 5%. 

In order to directly measure the r{t) curve, the EMuS sys- 
tem was run on the surface at the MIT-Bates Linear Acceler- 
ator Center in Middleton, MA. Candidate muon tracks are se- 
lected if they pass through two adjacent chambers on two paral- 
lel planes. A series of data cleaning cuts are applied to remove 
hit pairs created by noise and accidental triggers. Since the po- 
sitions of the wire chambers that fire are known, an estimate of 
the angle of the muon trajectory (8) can be calculated. Once the 

SNO X-Y Coordinate 0.53 cm 

SNO Z Coordinate 0.32 cm 

Floor Level* 0.17 cm 

Wire Placement 0.08 cm 

Wire Spacing* 0.18 cm 

Gaps Between Modules* 0.14 cm 

Time to Radius Conversion 0.28 cm 



Overall 



0.74 cm 



angle is known, the radii of closest approach are related as: 

R\+R 2 =Dcos8 (6) 

where D is the distance between each wire. A trial r(f) function 
(p(f) = at 2 + b) is used to estimate R 2 as a function of the time 
from the other chamber: 



R' 2 = Dcos0-p(*i) 
A least-squared parameter B is constructed 
B = (p(t 2 ) - R' 2 ) 2 



(7) 



(8) 



Table 2: Uncertainties associated with wire positioning. Uncertainties marked 
by an * do not apply to all wires. 



and then minimized. The resulting r(t) curve is shown in Fig- 
ure [7] Slices in time show a Gaussian shape, where the maxi- 
mum width of these slices is 0.24 cm, which is taken as the un- 
certainty on the time-to-radius conversion. The fit also extracts 
a negative time offset of 70 ns, which is caused by delays intro- 
duced by the electronic signal chain. This time offset slightly 
decreases the efficiency for reconstructing events, but does not 
significantly change the reconstruction accuracy. Running con- 
ditions varied slightly between Bates lab and underground at 
SNO (mainly due to ambient pressure and operating voltage) 
and simulations were used to correct for these changes. The 
extrapolation provides an additional uncertainty of +0.14 cm, 
yielding a total uncertainty of +0.28 cm on the time-to-radius 
conversion model. 

5. Data Selection 

A number of data quality checks were made to find candidate 
muons that went through both SNO and the EMuS system. Six 
of the EMuS wires were removed from the analysis because of 
their abnormally low or high trigger rates. A small number of 
channels had multiple recorded hits in a single event. For such 
events, only the first hit in time was considered part of the muon 
track reconstruction algorithm. 
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EMuS event level cuts were defined to select muon events 
throughout the run of the experiment. A minimum of three wire 
planes had to fire in order to ensure proper reconstruction. The 
event also had to have fewer than 30 wires fired so as to re- 
duce contamination from electrical pickup. Finally, runs with 
increased human activity above the detector, due to calibrations 
or source manipulation runs, were removed from the data anal- 
ysis. A total of 62 EMuS events passed all run selection criteria. 

To correlate these candidate events with the SNO detector, all 
of the relevant SNO runs were examined with an event viewer. 
Of the 62 EMuS events, 32 corresponded to a muon track pass- 
ing within the volume of the detector confined by the PS UP 
structure, while 16 corresponded to an event where a muon 
passed external to SNO's PMT support structure and was there- 
fore seen only by the outward looking PMT tubes. The remain- 
ing 14 EMuS events did not traverse the cavity. Of the 32 muon 
tracks within the SNO detector volume, 30 were properly re- 
constructed by SNO's muon fitter. The EMuS system ran for 
94.6 days of livetime, giving a rate of 0.32 reconstructed coin- 
cident events per day. 

6. EMUs Reconstruction 

By utilizing tracks that reconstruct in both SNO and the 
EMuS system, one can determine the final muon track recon- 
struction accuracy. A Monte Carlo-based method is used to de- 
termine such reconstruction characteristics. For each real data 
event that is reconstructed in both the SNO and EMuS detector, 
a series of random test tracks are generated. These Monte Carlo 
generated random tracks use the muon track as reconstructed by 
the SNO detector alone as a seed track, but its vertex and direc- 
tion are allowed to vary; with up to 66 < 10° variations in re- 
construction angle and up to 6b^ < 100 cm variations in impact 
parameter. Subsequently, these generated Monte Carlo tracks 
are then compared to the hit pattern as recorded in the EMuS 
tracking chamber. The negative log likelihood value (hereafter 
referred to as the likelihood) for each generated track is calcu- 
lated to determine the overall compatibility of the SNO muon 
reconstruction algorithm with tracks reconstructed in the EMuS 
system. The likelihood is given by the following functional 
form: 




20 40 

Impact Parameter Difference (cm) 



Figure 8: Angular difference vs impact parameter difference between SNO's 
muon fitter and the EMuS system for one event. The color scale indicates the 
density of possible tracks weighted by their likelihood. 



parameter difference) without significantly altering the hit pat- 
tern recorded by the EMuS system. 

Since this ambiguity exists only in the EMuS system and not 
in SNO's muon tracking algorithm, we can compare tracks re- 
constructed in the two systems by assuming either (a) the im- 
pact parameter is fixed or (b) the reconstructed track direction 
is fixed. To test the validity of these assumptions, an ensemble 
of fake data sets is generated both with and without accounting 
for track correlations in the EMuS system. The results from 
these Monte Carlo tests are shown in Figure [9] Correlations 
have no effect on the angular mis -reconstruction or the means 
of the distributions, but they do broaden the impact parameter 
mis-reconstruction by as much as 10 cm. We conclude that the 
EMuS-SNO tracks are sensitive enough to constrain the angu- 
lar reconstruction and impact parameter bias of the SNO muon 
fitting algorithm, but not the resolution of the impact parameter 
reconstruction. 



[bi-p(tdf 



(9) 



Figure 10 shows the results of applying the two assumptions 
to the 30 reconstructed EMuS-SNO events. The data are fit- 
ted to the functional forms of Equations [4] and [5] Due to the 
small number of events, the weights and relative widths of the 
secondary gaussians are fixed to their values from the earlier 
simulations. We find that the angular width is 0.61° + 0.06°. 
The impact parameter bias is 4.2 ± 3.7 cm, while fit impact pa- 
rameter width is 18 + 11 cm. 



7. Conclusions 



where b; is the impact parameter between the simulated track 
and the i th wire, p(f,) is the expected radius given the TDC time 
recorded for the wire and cr, is the wire position uncertainty. 
Wire hits that reconstruct at greater than 5cr from the main track 
are essentially removed to avoid reconstruction bias. 

Figure|8]shows the most likely tracks for a single event based 
on this method. The distribution is the projection of a cone, and 
indicates that there is a degeneracy between the angle and track 
reconstructed by the EMuS system. This is expected because 
if the track direction is changed (raising the angular difference) 
the placement of the track can be changed (raising the impact 



The combined data from the SNO detector and the Exter- 
nal Muon System have demonstrated that the SNO muon re- 
construction algorithm is accurate to the level needed by the 
neutrino-induced atmospheric flux analysis. The EMuS anal- 
ysis places a constraint on the angular reconstruction to better 
than 0.61° + 0.06° and on the impact parameter bias to better 
than 4.2 ±3.7 cm. The latter constraint is in good agreement 
with other methods using cosmic -ray data in SNO JT). We be- 
lieve the method employed here is a unique, low-cost way to 
explicitly verify the validity of muon track reconstruction for 
deep underground experiments. 
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Two Gaussian Angle Fit Without Simulating Correlations 
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Two Gaussian Impact Fit Without Simulating Correlations 
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Figure 9: Results from fitting the angular (left) and impact parameter (right) distributions of the ensemble of the generated simulated data sets according to 
Equations[4]and[5] respectively. The top plots show the results of fitting the distributions directly from SNOMAN Monte Carlo simulation package without taking 
into account correlations between angle and impact parameter reconstruction in the EMuS data. The bottom plots show the results with the inclusion of these 
correlations. 
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Figure 10: Gaussian fit to the data jointly reconstructed by the EMuS-SNO systems. Figure shows both angular (left) and impact parameter (right) difference. 
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